Failure of a granular step 
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We investigate the gravity driven rapid failure of a granular step composed of non-cohesive steel 
beads. The step is initially held together with electromagnets, and released when the current is 
switched off. We visualize the surface and the motion of the grains during the entire relaxation. 
The initial failure occurs at the surface and the subsequent flow is also confined to the surface as the 
step relaxes to its final state. The final shape of the surface is almost linear, depends on the initial 
angle of the step, and is not sensitive to the size of the grains. The average final slope of the pile is 
only slightly lower than the angle of repose of a pile formed by slowly pouring particles on to a flat 
surface. The evolution of the step is compared with a proposed convective-diffusion model of our 
system. The qualitative features of the relaxation are captured by the model after a flow-dependent 
dissipation parameter is introduced. 

PACS numbers: 45.70.-n, 45.70.Ht, 81.05.Rm 
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I. INTRODUCTION 

How a granular pile fails is important in understand- 
ing a diverse range of physical phenomena such as land- 
slides, breaking of retaining walls and even sand cas- 
tles Q, 0, ET- The nature of the failure depends on 
the applied constraints and the rate of failure. For ex- 
ample, in case of a retaining wall which holds granular 
matter, Coulomb's method of wedges is used to analyze 
the failure when the wall shifts slowly and the material 
undergoes quasi-static failure Q- This analysis tech- 
nique assumes a planar slip surface inclined at an angle 
which is determined by the weight of the material, fric- 
tion between the grains and the confining force exerted 
by the retaining wall. If the failure is catastrophic and 
the wall completely disintegrates, the materials acceler- 
ates rapidly and the resulting evolution can be rather 
different. 

Modeling rapid failure is challenging because of the ex- 
istence of a broad range of contact times between parti- 
cles, the transient nature of the dynamics, and a lack of a 
robust granular theory which encompasses both dynamic 
and static regimes. The situation is further complicated 
by a lack of quantitative data. A full scale simulation 
is difficult because of the wide range of grain-grain con- 
tacts that can exist during the failure. In the last decade, 
a continuum model has been developed by Bouchaud, 
Cates, Ramprakash, and Edwards (BCRE) for thin sur- 
face flows p . This approach greatly simplifies the com- 
plexity of the problem by assuming a static region and 
a thin dynamic region flowing with constant speed, and 
yields a convective-diffusion equation which can be then 
solved with appropriate boundary conditions. The model 
has excited much theoretical interest 0, S @i an d nas 
been extended to thick flows as well 0. 

Boutreux and de Gennes [|| have applied the BCRE 
equations to the evolution of a granular step. They pre- 
dict a final angle of inclination of the surface which is 
lower than the angle of repose for a pile formed by pour- 
ing grains on a flat surface with a low rate. However, few 



observations of internal flow exist, and the application of 
the BCRE approach to this situation is based largely on 
conjecture. 

Recent experimental investigations have studied such 
a situation when a cylindrical pile is released as a re- 
taining wall is rapidly lifted 0, . These studies have 
focused on the dependence of the shape of the surface on 
the preparation of the pile, and the aspect ratio of the 
pile. However, the flow itself has not been well character- 
ized which is important in developing an understanding 
of rapid failure. For instance, one would like to know 
where the initial failure occurs, i.e. does it develop at 
the surface or at depth. Furthermore, knowledge of the 
velocity fields inside the pile is important to determining 
the dissipation in the system. 

In this paper, we report an experimental investigation 
of the relaxation of a granular step using high speed imag- 
ing to address some of these open issues. Our technique 
allows us to measure the motion of the grains inside the 
pile. A key innovation is the use of magnetic fields to 
hold the initial step in place which results in a well con- 
trolled release with zero initial velocity. We find that the 
flow during the relaxation is confined mostly near the 
free surface. Our experiments show that the final angle 
of repose does not depend strongly on the initial angle 
of inclination, and in fact is similar to a pile obtained by 
pouring particles at a low rate. However, the detailed 
surface profile of the surface varies slightly for each case. 
For a pile formed by pouring at a low rate, we observe 
a linear surface. But for a step, the final profile shows 
a logarithmic deviation from linearity. We also compare 
the data with the model of Boutreux and de Gennes 
and find qualitative agreement after introducing a flow- 
dependent dissipation parameter. 



II. EXPERIMENTAL APPARATUS 

Figure m a) shows an example of a granular step. The 
step is created with steel beads between two parallel 
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FIG. 1: (a-h) A sequence of velocity fields of the grains su- 
perposed on images of the step after its release at t = Os. 
Velocity fields are measured by direct particle tracking, (a) 
t = 0.0 s, (b) t = 0.04 s, (c) t = 0.085 s, (d) t = 0.125 s, and 
(e) t = 0.175 s, (f) t = 0.25 s, (g) t = 0.45 s, (h) t = 1.50 s. 
The direction and length of the vectors corresponds to the 
mean velocity of the particle at that location (see text). 
(d = 2.8 ± 0.2 mm, h = 15cm). 

plates separated by width W = 2.5 cm. The pile is ini- 
tially held in place with the force generated by a set of 
electromagnets placed behind the apparatus. The top of 
the pile is flat and the right side of the step has an ini- 
tial angle 8i = 90 degrees. The height of the pile h is 
15 cm, in this case. The length of the step is such that 
the step can be considered to be semi- infinite in extent. 
Beads with diameters d = 0.9 ± 0.1 mm, 1.7 ± 0.2 mm, 
and 2.8 ± 0.3 mm that are spherical to within 20% are 
used in the experiments. A distributed filling procedure 
is adopted to uniformly build the pile. The particles are 
randomly packed and do not show crystalline order even 
near the side walls due to the irregularity in the shape 
and the size distribution of the particle. 

The relaxation of the step is actuated by switching off 
the electromagnets, and we have checked that the par- 
ticles do not show any effects of residual magnetization. 
This technique allows us to release all the particles simul- 
taneously. Furthermore, the particles are not susceptible 
to shearing that can occur while using a mechanical gate. 
A layer of 2.8 mm beads is glued to the bottom surface 
over which particles move as the pile relaxes. We image 
the relaxation of the pile through the transparent glass 
front wall using a mega-pixel resolution digital camera 
with a frame rate of 2,000 per second. The entire pile 
is imaged with sufficient resolution to track the motion 
of individual particles to sub-pixel accuracy using a cen- 
troid technique 0. Thus we obtain the shape of the 
surface and the velocity field during the entire relaxation 
process. 



III. THE LOCATION OF FAILURE 

The failure of a step is shown in Fig. ^ The initial pile 
is observed to deform rapidly initially, and then comes 
to rest with the free surface inclined at an angle to the 
horizontal. It can be noted that the surface is clearly non- 
linear at early times, with the top right corner remaining 
recognizable as it falls through more than half the total 
height. After this time, the surface appears more or less 
linear but systematic deviations can be also noted. 

Superposed on the shape of the pile are the mean ve- 
locity fields of the grains at that instant. Here the data 
is plotted for clarity after measuring the grain motion 
over a 0.5 ms time interval and averaging the obtained 
velocities over a 2. Id by 2. Id grid. It can be seen that the 
grains near the free surface flow with the greatest veloc- 



FIG. 2: (a) The velocity V as a function of depth z for 0i = 90 
degrees. The dashed lines shows the increasing phase with 
t = 0.01 s, t = 0.05 s, t = 0.085 s, and t = 0.125 s and the 
solid lines shows the decreasing phase with t = 0.16 s, t — 0.2 
s, t = 0.29 s, and t — 0.44 s. The velocity is observed to 
decrease smoothly to zero, (b) The maximum velocity at the 
surface v s as a function of time. The dashed line has a slope 
equal to the acceleration due to gravity, (h = 14 cm), (c) 
The corresponding depth of the flow. 



FIG. 3: (a) A plot of the final shape of the surface as a func- 
tion of the horizontal axis x and the vertical axis y normalized 
by the height of the pile. For a poured heap, we observe a 
linear profile, whereas the shape for the step with 9i — 90 de- 
grees is described by the equation ay/h + bln(y/h) +c = x/h. 
Inset: the angle of repose of the pile as a function of the width 
of the container, (b) The horizontal displacement versus the 
vertical displacement of the center of mass of the pile for var- 
ious combinations of step height and particle diameters. 



ity. In fact, soon after the magnetic constraint is removed 
at t = s, it can be noted that the velocity decreases to 
zero within a few grain diameters as one moves from the 
surface to the bulk [see Fig.QJb)]. Thus at early times, 
the flow is localized near the free surface and the location 
of failure is very different from that for a quasi-static re- 
moval of the confining wall. At later time the flow gets 
deeper but is still confined mostly to the surface. 

To clarify the magnitudes of the relative velocities as a 
function of distance from the surface, we plot the mean 
velocity V recorded at the center of the surface (h/2) in 
Fig-Efa) as a function of depth z normal to the surface at 
various time instances. At each time instance, the veloc- 
ity appears to decay approximately linearly with depth. 
At the earliest times, barely any motion is observed for 
depths greater than a couple of grain diameters. As time 
progresses and the step relaxes, the overall speeds in- 
crease rapidly before reaching a maximum and then de- 
crease somewhat slowly back to zero. 

The corresponding speed of the beads at the surface 
and the depth at which the velocity decreases to zero is 
plotted as a function of time in Fig. |2fb) and Fig. Htc), 
respectively. The data corresponds to a point approxi- 
mately mid- way up the inclined surface. A line with a 
slope corresponding to that for gravitational acceleration 
is also plotted. It can be noted that for 0j = 90 degrees, 
the initial acceleration of the beads is similar to gravi- 
tational acceleration but decreases before turning nega- 
tive as dissipation becomes increasingly important. From 
these plots, we also observe that the depth of flow grows 
initially reaches a maximum and then decreases. We have 
also added data corresponding to 6i = 45 degrees to test 
the effect on the initial angle on the relaxation. The 
overall behavior is similar except that the velocities are 
lower, and the peak occurs at a later time because the 
gravitational acceleration is lower. 
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IV. EVOLUTION OF THE FREE SURFACE 

We first discuss the final shape of the pile before ex- 
amining the change in its shape during the relaxation. 
The final surface profiles are plotted in Fig. E2a) for 
three different initial conditions. The profiles obtained 
for 6i = 90 degrees, and 0i = 45 degrees are plotted along 
with the shape for a pile formed by pouring the grains 
with a slow rate of about 20 grains per second inside 
the same container in which the relaxation experiments 
are conducted. The grains are poured from a point only 
slightly above the surface pile to prevent the particles 
from gaining substantial kinetic energy before they im- 
pinge on the pile surface. Because the surface flow can be 
intermittent when the grains arc poured slowly, we exam- 
ine the surface profile for the pile right after an avalanche 
event has occured to be consistent. The purpose of these 
second kind of experiments is to provide a baseline to 
compare the effects of the rapid motion, and to estimate 
the effects of the side walls on the observed shapes. As 
can be noted from Fig.[SJa), the steps relax only slightly 
further than the piles formed by slow pouring. The re- 
laxation for — 90 degrees is systematically greater than 
for 6 = 45 degrees. 

We find that the pile surface formed by slow pouring 
is described by a straight line. By contrast, the final 
shape after a rapid relaxation from an initial step shape 
is best described by an equation containing a linear and 
a logarithmic term: x/h = ay /h + 6 In (y/ti) + c, where 
y is the height of the surface at a horizontal position 
x, and a, 6, and c are constants. For the case shown 
in Fig. Eta), a = -1.865, b = -0.050, and c = 1.877, 
for 0i = 90 degrees, and a = -1.764, b = -0.014, and 
c = 1.779, for 6i — 45 degrees. We also varied the height 
h of the pile between 8 cm and 20 cm and found that 
a, 6, and c remains unchanged with standard deviations 
within 0.075, 0.021, 0.083, respectively. 

Our results for the shapes of the piles formed by slow 
pouring are inconsistent with the conditions reported in 
previous reports [l2L Il3j . In these previous reports, log- 
arithmic deviations from linearity were reported for piles 
formed by pouring and compared with models where dy- 
namics was neglected. However, we observe logarithmic 
deviations either after a rapid step relaxation or when 
grains are dropped on to the pile from a height so that 
they have substantial kinetic energy when they impact 
the pile. Thus we observe logarithmic deviations only 
when particles have substantial kinetic energy during the 
relaxation. 

In order to examine the effect of the side walls, we 
measured the surface of piles formed by pouring beads 
inside containers with various combinations of widths W 
and bead sizes. The surface was observed to be linear in 
all cases. The corresponding angle of repose as a function 
of the width of the system normalized by the particle 
diameters is shown in the inset to Fig.^a). The dashed 
line is an exponential fit to the data. The asymptotic 
value for the angle of repose is 25.5 degrees, and thus 



side walls do have a small effect on the measured surface 
profiles, but do not alter the qualitative features. 

As noted earlier, the final shape of the pile does not 
appear to depend significantly on step height. To illus- 
trate the scaling of the relaxation, we plot the horizontal 
Ax versus vertical Ay displacement of the center of mass 
of the pile with data obtained for three different parti- 
cle sizes [see Fig.^b)]. We plot the displacement of the 
center of mass rather than the length of the relaxation 
because it is more suseptible to fluctuations in the tail of 
the relaxation. A linear dependence is observed. Thus, 
it appears that the surface profile is not very sensitive 
to the grain size, which gives support to a continuum 
approach to describing the relaxation. 



V. COMPARISON WITH A 
CONVECTION-DIFFUSION MODEL 

As first discussed in the introduction, a phenomeno- 
logical continuum description of this situation has been 
developed by Boutreux and de Gennes [|| based on the 
BCRE model. They assumed that the flow during re- 
laxation is confined to the surface consistent with our 
observations. We next briefly reproduce their argument 
for the rate of change of the surface as the step relaxes 
before comparing the predictions with our data. They 
first assumed that the surface is approximately linear, 
and given by y = xtan9. Then it was argued that the 
rate of change of the total height of the pile is given by 



dt 



dR 

dx ' 



(1) 



where, R is the height of the liquid-like mobile phase, and 
v is the mean velocity of particles in the mobile phase. 
Next they assumed that R has a maximum Rq at the 
center and decreases smoothly to zero at the front and 
back edge of the step and therefore 



dR 

dx 



a- 



Ro tan 9 



(2) 



where, a is assumed to be a non-dimensional number 
with an absolute value of order 1 near the extremum and 
less than 1 elsewhere. Next they assume that the rate 
of change of the mobile phase at the maximum mobile 
depth Rq at the step-center is given by 



dR 
~dT 



7(tan0 — tan6>d)i?o, 



(3) 



where, 7 is a frequency which is related to the flow veloc- 
ity v and the diameter of the particle. They assumed 8d 
to be the angle of repose of the particles 6 r . This equa- 
tion comes from the BCRE model 0, and it is further 
assumed that convection has negligible effect on the rate 
of change of the rolling phase at its maximum. 

Next, they substituted Eq. J5J) in Eq. QJ to get R 
which in turn is replaced in Eq.10 to give a second or- 
der differential equation for the distance traveled s after 
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FIG. 4: (a) The evolution of the average inclination angle of 
the granular free surface, (b) The normalized runoff distance 
of the step cot 8 = s/h as a function of time. The dashed lines 
are the numerical integrations of Eq. (4) (see text), (c) The 
rate of change of the normalized runoff distance cot 9 = s/h. 
The dashed lines are fits to Eq. (4). The maximum occurs 
at higher angles for the higher initial angles, and seems to 
be independent of the initial height of the step, (d) Velocity 
dependence of the effective dissipation for the particles at the 
surface. 

rewriting the variable y in terms of s. This equation 
upon integration yields the rate of change of the runoff 
distance s normalized by the initial height of the step 
cot 9 = s/h 0: 

1- rfcot =ln(cotfl) -tanfl d cotfl + fc, (4) 
7 at 

where, k is a constant of integration. Equation Q pre- 
dicts a maximum for at 9 = 9d which can be tested 
against our experimental results. 

Figure0{a) shows the plot of the average surface angle 
of the pile corresponding to a step with three different 
initial angles 9i — 45 degrees, 60 degrees, and 90 degrees. 
Because the surface is not quite linear as discussed ear- 
lier, we average the angle over the entire free surface but 
exclude 3 particle diameters from the top and bottom of 
the step. The data shows a rapid decreases at early times 
with a crossover to a slower decrease to the final angle. 
In Fig. Efb) , we plot the corresponding normalized runoff 
distances cot 9 as a function of time. Then to compare 
the data to the model, the rate of change of cot 9 for the 
same data is plotted in Fig. 0Jc). The peak occurs well 
before 9 r is reached in all cases and depends on the initial 
angle. Thus if 9d = 9 r as assumed in the discussion of 
the model above, then even qualitative agreement is not 
seen. 

Next we allow 9d to vary and find that Eq. 0} can 
indeed qualitatively describe the data. In Fig.01c), the 
dashed lines are the fits to Eq. ©. The corresponding 
fitting parameters are listed in Table [fl This possibility 
was also discussed by Boutreaux, et al. while extend- 
ing the BCRE equations to thick flows. They anticipated 
that 9 r in Eq.[3]may have to be replaced by a higher neu- 
tral angle. Our observation of an higher 9d is also con- 
sistent with experiments on heap flows where an increase 
in neutral angle is observed with flow rate • 

We note from Table [I] that 9d increases systematically 
with 9i which results in higher flow speeds as can be 
noted from Fig. Hfb). Now we can also use the fitting 
constants to numerically integrate Eq. and obtain the 
normalized runoff distance as a function of time. The cor- 
responding curves (dashed lines) are shown in Fig. QJb) 
and qualitatively captures the evolution of the normal- 
ized runoff as well. Repeating the experiment with dif- 
ferent heights for the same 9i gives similar 9d- Thus we 
further confirm the scaling behavior of the system. 



6, (deg.) h (cm) 7 (s- 1 ) 6 d (deg.) k 



90 


17.5 


1.89 


68.4 


3.947 


60 


17.5 


15.26 


43.8 


1.156 


45 


12.5 


50.48 


36.2 


0.735 


90 


14.0 


3.92 


67.6 


3.820 


60 


13.0 


19.72 


43.9 


1.154 



TABLE I: Parameters 7, 6d, and k obtained by fitting Eq. gj 
to the data. The first three rows corresponds to the data 
plotted in Fig. EJc) . 



The fact that 9d is greater than 9 r and increases with 9i 
seems to be self-consistent with the idea that dissipation 
in granular flows occurs not only because of sliding fric- 
tion but also because of inelastic collisions. Figure ^d) 
shows the velocity dependence of the normalized effec- 
tive dissipation force at the surface layer. This has been 
obtained by balancing the average forces on a particle 
due to gravitation acceleration along the inclined surface 
and the dissipative force normalized by the weight of the 
material fdis at the surface of the pile and can be simply 
calculated from fdis = sin. 9 — a/g. One notes that the 
dissipative force is an increasing function of velocity and 
has a maximum at the maximum velocity before finally 
dropping down to the static value as the material comes 
to rest. The observed hysteresis arises possibly because 
of differences in the contribution of sliding friction during 
the deceleration and acceleration phase due to differences 
in inclination angles. However, without a detailed dissi- 
pation model for transient flows it is difficult to interpret 
the data further. 



VI. CONCLUSIONS 

In summary, by visualizing the failure of a granular 
step, we have found that it relaxes to angles only slightly 
less than the angle of repose of a pile formed by slow pour- 
ing. Furthermore, the flow is observed to be confined to 
the free surface with the velocity decaying smoothly with 
depth. The step is observed to relax rapidly at first with 
the maximum rate occuring at an angle greater than the 
angle of repose. Thus the resistance to motion offered 
by inelastic collisions in addition to sliding friction must 
be important in determining the evolution of the pile. 
From our experimental observations, it appears that a 
convectivc-diffusion model can describe the qualitative 
features of a rapid step relaxation by introducing a ve- 
locity dependent dissipation parameter. It is interesting 
that this is the case considering that some of the as- 
sumptions in the model such as constant flow speed are 
not met in the experiments. Nonetheless, the fact that 
such an approach can capture the dynamic evolution is 
interesting considering the difficulty in modelling dense 
granular flows. Further work is required to develop and 
understand the physical implication of the model now 
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that experimental data is available. 

Note added - During the review of our manuscript we 
became aware of two papers which discuss granular col- 
umn collapse in two and three dimensions [la . Il6l | in ad- 
dition to Ref 10. These studies examine the scaling of the 
runoff distance with column width and may be considered 
complementary to our study focusing on semi-infinite ex- 
ten columns. 
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